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ABSTRACT: Hamiltonian and Schrodinger evolution equations on finite- 
dimensional projective space are analyzed in detail. Hartree-Fock (HF) manifold 
is introduced as a submanifold of many electron projective space of states. Evo- 
lution equations, exact and linearized, on this manifold are studied. Comparison 
of matrices of linearized Schrodinger equations on many electron projective space 
and on the corresponding HF manifold reveals the appearance in the HF case a 
constraining matrix that includes matrix elements of many-electron Hamiltonian 
between HF state and double excited determinants. Character of dependence of 
transition energies on the matrix elements of constraining matrix is established by 
means of perturbation analysis. It is demonstrated that success of time-dependent 
HF theory in calculation of transition energies is mainly due to the wrong be- 
havior of these energies as functions of matrix elements of constraining matrix in 
comparison with the exact energies. 
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Introduction 

There exist two widely used by quantum chemists simple approaches for 
calculation of excitation energies of many electron systems having at their 
heart the Hartree-Fock (HF) theory. In both of these approaches it is presup- 
posed that optimal HF molecular spin orbitals (MSOs) are already calculated 
and that the excited many electron wave functions are linear combinations 
of determinants obtained from the HF one by all possible single excitations. 
The first, most simple, approach is based on the so-called single-excitation 
configuration interaction (CIS) method and can be technically described as a 
diagonalization of projection of many electron Hamiltonian on the subspace 
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of single excited determinants. In the second approach, where the so-called 
time-dependent (TD) method is used, the same subspace of single excited de- 
terminants appears as a tangent space to the HF manifold. However, arising 
on this subspace operator does not coincide with the Hamiltonian projection 
but involves parameters, accounting indirectly for the double excitations from 
the HF state. 

TDHF equations were derived by Dirac from his time-dependent varia- 
tional principle by constraining trial wave function to remain a single de- 
terminant at all times Interpretation of linearized TDHF equations in 
terms of harmonic oscillations of a certain fictitious many particle system in 
a neighborhood of the HF minimum was described by Thouless [2]. More 
rigorous analysis of TD theories based on methods of modern geometry was 
undertaken by Rowe and his co-authors [3]- [6]. He suggested to treat con- 
strained time- dependent many-body quantum mechanics as a Hamiltonian 
system on a symplectic manifold. This system can be obtained from Dirac's 
extremal condition of an action integral. 

Our approach, presented here, is closed in its concept to that of Rowe. 
We, however, simplify analysis of TD theories by fixing from the very begin- 
ning a concrete altlas covering symplectic manifold M and work within one 
selected chart of this atlas centered at the HF optimal determinant. This 
strategy allows us to use, after realification, the symplectic structure of the 
parameter space instead of usually more complicated symplectic structure of 
a manifold itself. Then we write down the first order differential equation 
on a symplectic parameter space. The right-hand side of this equation is 
supposed to be a symplectic gradient vector field corresponding to a local 
representative of a smooth real-valued function on M. Any such function is 
called Hamiltonian function in commonly accepted in theory of symplectic 
manifolds terminology. The corresponding differential equation is also called 
Hamiltonian. Subsequent complexification of the realified parameter space 
leads to the Schrodinger-type evolution equations. 

To make our reasoning as independent as possible from numerous text- 
books on modern geometry, in Section II necessary definition from manifold 
theory together with few simple examples are given. Section III is dedicated 
to TD theory on finite-dimensional projective spaces. The corresponding 
linearized Hamiltonian and Schrodinger equations should be considered as 
finite-dimensional versions of the exact evolution equations. In Section IV 
HF manifold is introduced and its several atlases are described. Linearized 
and exact evolution equations on HF manifold for arbitrary quadratic Hamil- 
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tonian function are derived. Then the analogous theory is developed for 
arbitrary (not necessarily quadratic) energy functions depending on matrix 
elements of 1-density idempotent operators (Hamiltonian functions on Grass- 
mann manifolds). 

Transition energies obtained with the aid of TD theories are usually more 
close to the exact ones than, say, the CIS energies. In Appendix A it is 
demonstrated that this is an effect (somewhat paradoxical) of a wrong be- 
havior of TD transitions energies as functions of certain complex parameters 
in comparison with the behavior of the exact energies. In Appendix B ex- 
pressions for derivatives of the Gram-Schmidt parametrization function are 
collected. To the best of our knowledge, this parametrization was introduced 
in quantum chemistry by Garton [7]. 



Basic Definitions 



Many methods of quantum chemistry are based on a simple idea of 
parametrization of a certain subset M of state vectors from the p— electron 
sector of the Fock space or from the corresponding projective space by el- 
ements of some parameter space of a relatively small dimension with sub- 
sequent optimization of chosen parameters using one or other optimality 
criterion. A set M may be a surface or, more generally, a manifold or a 
variety. We will give all definitions and discuss most important properties 
of some typical set M supposing that it is either locally Euclidean or locally 
Hermitian space, which means that each its point has a neighborhood home- 
omorphic to an open subset of the number space K" where K is either the 
field M of real numbers or the field C of complex numbers. 

Let us start with general definitions not presupposing at first, that the 
set M is embedded in some space of states. All this definitions may be found 
in many excellent books dedicated to the manifold theory (see, e.g., [5]-[TB]) 
including books written by mathematicians especially for physicists [T^ [TB] . 

Definition 1. Triple c = {U, (f, K"^) is called a chart on M of dimension m 
if (1) U is a subset o/M; (2) (f is a bijection of U on an open set in K™. 

Subset U is called the domain of the chart c and mapping is a (local) 
coordinate system on M. The inverse of if is called a local parametrization 
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of M. If X E U then c is said to be a chart on M at point x. If ip{x) = then 
c is called a chart centered at x. 

Let c = {U, (fi, K"") and c' = ([/', (^', K™) be two charts on M. These charts 
are called compatible if (1) both ip{U fl U') and (p'{U H U') are open in ffiC^; 
(2) the mappings (that are called transition functions) ^p' oip~^ : ip{U fif/') — *• 
^'{U n U') and o (^'-^ : ^'{U H f/') ^ <^(f/ n f/') are K-analytic. 

Definition 2. Family A of charts {(f/j, y^j, K™)}j called an analytic atlas 
o/M z/(^ijM = |J [/j/ f^Sj 74ny two charts from A are compatible. 

i 

Definition 3. The set M with an analytic atlas on it is called ^-analytic 
manifold. 

Definition [3] characterizes manifold with the aid of some its concrete atlas. 
More elegant definition may be given on the base on the notion of compatible 
atlases. 

Two atlases A and A! of M are called compatible if ^ U is an atlas of 
M. It is easy to show that compatibility of atlases is the equivalence relation 
on the set of all atlases of M (see, e.g.|9|). 

Definition 4. Class of equivalent atlases defines ^-analytic manifold struc- 
ture on M. 

If the ground field is M then the requirement of analyticity of transition 
functions is too restrictive. Instead the notion of differentiable manifolds of 
class C^{k = 0, 1, ... , oo) may be introduced. Transition functions of such 
manifolds are supposed to be continuously differentiable up to order k. It 
is clear that transition functions of class C'' are actually C'^-diffeomorhisms. 
Remind as well that R-analytic manifolds are C°° ones but not vice versa. 

On the same set different manifold structures may exist. The simplest 
standard example is the set M of real numbers: two charts (M, id, M) and 
(M, (f, M) where id : x ^ x and ip : x ^ x^ determine two different manifold 
structures on M. Indeed, each chart endows M with structure of C°°-manifold. 
But these two charts are not compatible: the transition function ip o id~^ : 
X ^ x^ is smooth and bijective but the inverse mapping id o ip~^ : x ^ xs is 
not differentiable at the origin. 

The notion of manifold is a very abstract generalization of the classic 
notions of opened smooth curves and surfaces in Euclidean spaces. If a : 
(a, b) —>■ M™ is a smooth curve, then it can be interpreted as a smooth 1- 
dimensional manifold with the chart c = (cr(a, 6), o"~^, M). Another example 
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is a graphic of smooth function / : R"^ M (the set of pairs (x, f{x))). It is 
m-dimensional manifold M/ in W^^^ with atlas consisting of a single chart 
(JVlf,{x,Xm+i) ^ x,W^). The corresponding parametrization mapping is 
X {xj{x)). 

Classic M-analytic manifold that can not be covered by a single chart 
is the unit sphere S™"^. Among its atlases probably the most simple is 
constituted by the charts Cae = {Uaey'^aey^"^^^) where e = ±1, 

U^, = {xe S""-^ : ex^ > 0} (1) 

and ipae '■ X ^ [xi, . . . ,Xa, ■ ■ ■ , Xm)- Here the hat over variable means that 
this variable is omitted. The inverse mapping (parametrization) is 

^"^^ ^ \5v/l-l|x(")p ifj=« 

where < 1. The minimal number of charts covering S™~^ is equal 

to 2. These charts may be constructed with the aid of, say, stereographic 
projection. 

Till now all manifolds we considered were actually subsets of Euclidean 
spaces. There arises a natural question: Do there exist manifolds that can 
not be realized as subsets of appropriate number spaces ? In the case 
of different iable manifolds the answer is negative. In 1936 Whitney proved 
that any different iable manifold of dimension m admits embedding into Eu- 
clidean space M^™"*"^ [17]. For C-analytic manifolds, however, the situation 
is completely different. It is easy to show that compact C-analytic manifold 
of positive dimension can not be embedded into Hermitian space C™. In- 
stead there exists a broad class of compact C-analytic manifolds that can be 
embedded into appropriate projective spaces. Some of such manifolds that 
occur in quantum chemistry, will be considered in detail in the next sections. 

Let M and K be two manifolds of dimension m and n, respectively, and 
/ be a mapping / : M — > 3\f. If (f/, (p, K"^) is a chart on M and (V, ^p, W) is a 
chart on such that f{U) C V then the mapping ipofoip'^ is a classic vector 
function of m variables defined on a certain open subset of K"^ and having 
subset of as its range. This mapping represents f in the charts under 
consideration. Mapping M K is called different iable at point x G M if its 
representative in selected charts {U, ip, W^) and (V, ip, K") is differentiable at 
point ip{x). It is easy to show that this notion of differentiability does not 



5 



depend on the choice of charts. Different iable on M functions can therefore be 
defined as functions differentiable at each point of M. Mapping f : M ^ "N 
is called morphism of manifolds if any its representative belongs to the class 

for C'^ manifolds and is analytic for analytic manifolds. 

Local linearization of a manifold in a neighborhood of some its point 
leads to very important notion of tangent space. Tangent spaces can be 
introduced in a number of equivalent ways of which we describe probably 
the simplest one (see, e.g., [H]). Let us consider the set of pairs (x, cr) where 
a; e M, cr : (— e, — >• M is differentiable mapping (curve) such that a{0) = x 
(curves passing through a point x G M). Two pairs (xi, cxi), (x2, o"2) are 
called equivalent if Xi = X2 and within some chart c = (f/, y?, K") on M at x 
the derivatives Do[v^°cri] and DQ[ipoa2\ coincide. Class [(x, cr)] of equivalent 
pairs is called tangent vector to M at point x and the set of all classes is 
called the tangent bundle of manifold M and is denoted TM. To put it 
more precisely, TM is a total space of the tangent bundle which is a triple 
(TM, TT, M) where M is a base and vr : [(a;, a)] ^ x is a projection of this 
bundle. Fiber t[~^{x) over point x G M is called tangent space to M at 
point X and denoted T^.M. It contains classes [(x, a)] with fixed x. Since 
each derivative -Do[v^ ° is, by definition, a linear mapping M — IC", we 
can introduce a mapping 6c ■ [{x, a)] Do[ip o a]{l) G that is obviously 
injective. Existence of differentiable curve a{t) = (f~^{(f{x) + tv) for any 
vector V G IK*" shows that 9c is actually a bijection. With the aid of vector 
structure from is transferred to the tangent space T^-M which becomes 
therefore m-dimensional vector space. Mapping 6c depends, of course, on 
the chart chosen. Indeed, if at point x another chart c' = (f/', y?', K™) is 
taken then 9c' = -D^(2^)[v?' o {p~^] o 9c where, by definition, D^p(^x)['^' ° '^~^] is 
an isomorphism of the vector space K™. This shows that the transferred 
vector structure on T^-M does not depend on the choice of concrete chart. 
In particular, vector space is a manifold that can be covered by a single 
('natural') chart (K™, irf, K™). As a representative of class [(x, a)] one can 
always choose a pair {x,x + tv) where v = Do[(p o cr](l). We have 9~^ : 
V [{x,x + tv)]. It is a common practice to identify TK™ with the direct 
product IK*" X and consider the tangent space T^W^ as the set of pairs 
(x, v) where f is a vector outgoing from point x (see, e.g., [13]). 

Let / : M — > !N be a morphism of manifolds, x G M and y = f{x) G K. 
Let us suppose that c = (f/, ip, K™) and c' = (V, ip, K") are charts on M at x 
and on K at y, respectively, and f{U) C V. The derivative of the mapping 
F = ip o f o (p~^ at the point (p{x) is a linear mapping K'" — > K". Using 



6 



the aforementioned bijections : T^M and O^' : Tj(a)!N" K" we can 

define linear transformation 6^' o D^(^^^F o 9^^ from T^-M to Tj(2,.)3\f which is 
called the tangent mapping to / at point x and is denoted as Txif)- This 
mapping does not depend on the choice of charts c and c'. In particular, any 
local parametrization 7 : K"* — > M of a certain open subsets of M induces 
tangent mapping T^i'j) '■ T^-K"* — > T^(j.)M at each point of the parameter 
space. 

In the case of morphism / : M ^ K" instead of the tangent mapping 
Ta;(/) the differential of / is used. Differential is a linear transformation 
dxf '■ Tj;M — > that is defined in the following way. If c is a natural chart 
on K" then d^f — Oc o Tx{f). It is pertinent to note, however, that the 
delicate difference between notions of tangent mapping to / and differential 
of / is usually ignored in mathematical literature and tangent mappings are 
also called differentials. We will also keep to this tradition. 

If / is a morphism M ^ K then differential dxf is an element of the 
vector space T^-M* dual to the tangent space to M at point x. This dual is 
called the cotangent space and its elements are the so-called co-vectors. 

It is well-known from linear algebra that there is no basis-independent 
(canonical) isomorphism between a vector space and its dual. An additional 
algebraic structure on vector space is required to perform canonical transfor- 
mation of vectors to co- vectors and back. Need in such transformation arises 
in almost all physical and many mathematical theories. Since transformation 
co-vector— vector will be used in subsequent sections, we found it reasonable 
to remind here the necessary definitions. 

Let g : E X E ^ Whe a, non-degenerate bilinear form on a real vector 
space E. It defines a canonical isomorphism of this space and its dual E*. 
Indeed, for any fixed v & E the partial mapping 

g{-,v) g{u,v) (3) 

is a linear functional on E and the mapping 

e^:v^g(;v) (4) 

is the aforementioned isomorphism. If {e^}i<^<n is a basis of E then the 
inverse mapping : E* ^ E may be written as 

n r n 

H=l \_v=l 
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where g^^ — g^^ is the matrix inverse to the matrix 

9i,v = 9{e^, e^) (6) 

of bihnear form g relative to the basis {e^}. 

When the ground number field is R, there are two most important par- 
ticular cases: bilinear form is non-degenerate symmetric or non-degenerate 
skew-symmetric. In the first case pair {E, g) is called Euclidean space, in the 
second case it is a symplectic space. Symplectic forms are usually denoted 
by the symbol ur. 

When the ground number field is C, functionals I : E ^ C such that 

l{u + v) = l{u) + l{v) and l{cu) = cl{u), u,v E E,c E C (7) 

are usually considered. The are called |-linear functional and the vector space 
(over C) of such functionals is also denoted E*. Instead of non-degenerate 
bilinear forms the so-called l|-linear (or Hermitian) forms are used. We 
consider Hermitian forms ^-linear with respect to the first argument: 

g{u, v) = ^ Uf^g^^v^ (8) 

IJ,,V=1 

This is consistent with widely used in physics Dirac's notations: {u\v) is 
fi-linear form |-linear with respect to u and {u\ is just ^-linear functional. 
Mathematicians usually prefer fi-linear forms that are ^-linear with respect 
to the second argument. 

Let M be a differentiable manifold of dimension m embedded in the Eu- 
clidean space M"-, / be a smooth mapping M" — >• M, and 7 : — > M be 
a parametrization of M, global if M is a surface and local if it is a man- 
ifold (without loss of generality one can always take M™ as a domain of 
parametrization mapping) . Stationary condition for the function / o 7 at a 
point X is 

dx [/ o 7] = di{.x)f o 47 = (9) 

Linear isomorphism dj;7 maps (parameter space) on the tangent space 
T^(j:)M which can be considered as a subspace of M". Image d^r/^ei) of the 
canonical basis vectors of the parameters space is a basis of the tangent space 
T!y(a;)M. Differential daf is a co- vector 

daf-J2Q^('')^y^ (10) 

j=i 
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where {dyj} is a basis of (M")* dual to the canonical basis {ej} of M". In 
more habitual for physicists Dirac's notations 

dJ = i2^(^)(^^\ (11) 

j=i 

For manifolds embedded in Euclidean spaces the differential of the mapping 
7 may also be conveniently written in Dirac's notation as 

d^^ = J2\£('^^)(^^\ (12) 

i=i * 

where {cj} is the canonical basis of the parameter space. 

Standard Euclidean scalar product g{ej,ej') = {ej\ejf) = 6jj' on may 
be used to identify this space and its dual. From Eq.Q it readily follows 
that 0~^((ej|) = \ej) and stationary condition Q) takes the form 

(e;i(rf^(,)/)M,7(e.)) = (2 = 1, . . . , m) (13) 

Geometrically this means that gradient of / (vector ©3^(^7(0;)/)) at a sta- 
tionary point 7(x) should be perpendicular to the tangent space T^M,. In 
quantum chemistry conditions of the type of Eq. f|T3|) are called Brillouin 
conditions. 

As has already been mentioned, tangent spaces to manifolds embedded in 
Euchdean spaces may be considered as a subspaces of the enveloping space 
and, consequently, they inherit its Euclidean structure (scalar product). In 
particular, the Gram matrix Gn'i^x) = {dx'j{ei)\dxj{ei')) (overlap matrix in 
the terminology accepted by quantum chemists) is defined on each tangent 
space, and, if smooth in x, endows M with the structure of Riemannian 
space. Riemannian metric is used to study the internal geometry of surfaces 
and manifolds. 

Let us return to our simple examples. Differential of the mapping 7:0;—*^ 
(x,/(x)) that parametrizes the graphic M/ C M"^"^^ of a smooth function 
/ : ^ M is 

dxj = ^\ei + -^{x)em+i){ei\ (14) 



=1 



Tangent space to M/ at point (x, /(x)) is spanned by the vectors dx'y{ei) = 
(^i + §^ii^)(^m+i and Gram matrix is Gw^x, /(x)) = 6ii> + §^^{x)^{x). Func- 
tion / may be written as pr^+i o 7 where prm+i{x,Xm+i) = Xm+i- Brillouin 
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conditions for this function (em+iMx7(ej)) = §^{x) = are just the classic 
stationary conditions for function / at point x. 

Differential of the parametrization mapping ([2]) is 



If / is a smooth mapping M."^ M then the stationary conditions for its 
restriction on the unit sphere S™~^ are 

= Ax, X G S"^-^ (16) 

dx 

m 

where ^(x) = ^ ^(x)ej is the gradient of /. In the case under consid- 

eration Brillouin conditions may be written in the form independent on the 
chart index. 

Let I be an open interval of M. It is a trivial manifold covered by a single 
chart {I,id,M.). Morphism cr M is called a smooth curve on M. Its 

differential dfCr is a linear mapping M To-(i)M. Such linear mapping is 
uniquely determined by the vector dtcr{l) which is called a tangent vector to 
curve a at point t E I. 

Vector field on M is defined as a mapping that to each point x G M puts 
into correspondence a vector from the tangent space T^M. For example, for 
any smooth function / on M the mapping x — >■ 0^(^)(c^x/) is a vector field 
on M where g{x) is smooth in x non-degenerate bilinear form on Ta;M. If ^ 
is a smooth vector field on M then the solution of the first order differential 
equation 

ait) = e(a(t)) (17) 

is called an integral curve of this vector field. 

In conclusion of this section it must be admitted that in applications 
the abstract manifold theory recedes in the background and the information 
given here is therefore somewhat excessive. The role of abstract theory re- 
duces to recognition of geometric object as a manifold and to selection of 
convenient local coordinates on this manifold. After the concrete atlas of 
the manifold under consideration is chosen and parametrization mappings 
are constructed, within a given chart instead of frequently complicated Rie- 
mannian or symplectic metrics on the tangent spaces one can successively 
use, as a rule much more simple, Euclidean or symplectic structures of the 
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parameter space. Examples of such a strategy are given in the next sections 
of the present work. 



Projective Spaces 



Projective spaces supply us with the simplest example of compact K- 
analytic manifolds (K = M or C). From physical viewpoint state of quantum 
system is a vector of the relevant Hilbert space determined up to an arbitrary 
phase pref actor being therefore a point of the corresponding projective space. 

The set of 1-dimensional subspaces (lines passing through the origin') of 
the vector space K""*"^ is denoted as Pn(lK) or as P(K"'"'"^) and is called the 
standard n-dimensional projective space over the ground field K. 

For any nonzero z € K""*"^ symbol \z\ stands for 1-dimensional subspace 
generated by vector z. Coordinates zi^^Z\^ . . . ^z^ of vector z are called ho- 
mogeneous coordinates of line \z\ (due to the property \z\ — \\z\ for any 
A G K\{0}). 

For each a = 0, 1, . . . , n let us define 

C/a = {N e P„(K) : ^ 0} (18) 

The mapping 

^^a '■ [-^0) ■ ■ ■ ) Z(x-\i ^ai -S^a+l) ■ ■ ■ i Zr\ ( , ■ ■ ■ i 1 — ■ ■ • ) ) G 

Zq, Zq, Zq, Zq, 

(19) 

is a local coordinate system on Ua and the family of charts Cq, = (?7q,, (pa, IK") 
is an atlas of K-analytic structure on P„(]K). The inverse mapping 

: (Co,...,Ca,---,Cn)^ [Co,---,l,---,Cn] (20) 

is a local parametrization of [/„ by elements of K". The hat over variable 
means that this variable should be omitted. 

There exists surjcctive mapping tt : ]K"+^\{0} Pn(K) defined by the 
relation t:[z) = [z]. Its restriction to the unit sphere 5*" for K = M and to 
the unit sphere 

S'-^' = {z e C"+^ : {z\z) = j2 ^^'^ = 1> (21) 

j=0 
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for K = C is also surjective. Since for any nonzero z G K""*"^ 

[2] = {+^,-^} if K = M (22) 

and 

\z\ n = {e''"^ : <^ G [0, 27r)} if K = C (23) 

it is possible to realize the projective space Pri(lK) either as a quotient of the 
unit sphere S*" modulo the equivalence relation 

z' ^ 7! = ±z (24) 

(K = M), or as a quotient of the unit sphere 5'^""*'^ modulo the equivalence 
relation 

z< ^z' = e'^z (25) 

(K = C). In further discussion we confine ourselves to the most interesting 
for us case K = C. It is easy to see that the mapping 

C = (Co, ■ • • , Ca, • • • , Cn) (Co, • • • , 1, • • • , Cn) 

1 I e'*^ 1 

"7rTlcy«°--'^--'^"'"y{^T^(Co.....^....c..)}(26, 

is a local parametrization of P„(C) realized as a quotient of S*^""*"^ modulo the 
equivalence relation fl25|) . Vector — ^===((^q, . . . , 1, . . . , ^„) is a representative 

of the corresponding equivalence class. In fact, we have M- analytic local 
parametrization 

7a : (Co, • • • , Ca, • • • , Cn.) ^ . • • • , 1, • • • , Cn) (27) 

V 1 + \CK) 

of the representatives of the equivalence classes (l25i) and 7q,(C") is a 2n- 
dimensional surface situated on S"^""^^. 

For geometric characterization of tangent spaces to the projective mani- 
fold Pn(C) it seems reasonable to start with the realification of the complex 
vector space C""^^ (see [T6]) that gives the real vector space CJ^"*"^ ~ R^"+^ 
with the standard basis 

Co, ei, 62, . . . , Cn, ieo, iei, ie2, • • • , ie^ (28) 
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Hermitian scalar product on C""^^ (which is supposed to be ^-hnear with 
respect to the first argument) may be written as 

{z\z') = {z,z')+\[z,z'] (29) 

where 

{z,z') = Re {z\z') (30) 
is the Euchdean scalar product and 

[z,z']=lm {z\z) (31) 

is the symplectic one, both of them are non-degenerate. Basis fl28l) is or- 
thonormal with respect to the Euclidean scalar product ( !30l) and is also stan- 
dard symplectic basis: [e^, Cj] = [iei,iej] = 0, and [ej,iej] = 5ij. Euclidean 
and symplectic scalar products are connected by the relation [z, z'\ = {iz, z'). 

Tangent vector space to 5^"^^ at some point z is the orthogonal com- 
plement to vector z in R^"'^^ with respect to the Euclidean scalar product 

T,52"+^ = {CgCS+^:(z,C)=0} (32) 

Point z G S'^"''"^ is a representative of the circle {e^'^'z : G [0, 2tt)}. Tangent 
vector to this circle at point z is ^e"^z|(^=o = i^- It is easy to ascertain that 

(^10 = ^ (^, C) = and {iz, C) = (33) 

As a result, the tangent space T[2]P„(C) at point [z] is isomorphic to (Cz)^ 
(the orthogonal complement to the line Cz with respect to the Hermitian 
scalar product on C""*"^). Orthogonal projection of arbitrary vector C, G C"^^ 
on the orthogonal complement to z is of the form 

If dzTT : T^C"^"'^ T[2]P„(C) is the differential of n at z, and r/i = d^iTl^i), 
V2 = dz'^{^2) ai'e two tangent vectors from T[2]P„(C) (their concrete nature 
is irrelevant) then it is possible to introduce on T[2]P„(C) the following Her- 
mitian scalar product 



Knimiz] = j-r-r^ (35) 







-(61 




16) 


(^1 


1^)2 
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Now let us consider parametrization (12 7p assuming that for each a the 
parameter space C" is embedded in C""*"^: C" C C" © Ccq,. Differential (i^7o, 
is an isomorphism T^M^" T^^(^)7q,(M^"). Simple calculations give 

= M /Lm^ - ^"(0)7.(0] (36) 

where ^ belongs to T^R^" ~ M^" and (7a (C); is the Euclidean scalar product 
(!30|) . It is easy to see that with respect to this scalar product vectors (136!) 
are orthogonal to 7a(C) and to ico,. Thus, with such an approach, 

= T,.(C)^'"'''®^7a(C) (37) 

and 

T7.(C)'^'"^' = T7.(C)7a(M2") © Mie„ (38) 

The expressions for the Euclidean and symplectic scalar products on the 
tangent space T^^(^)7q,(R^") are 

(a). . (ei,6)(i + iicf)-(ei,c)(c,6) ..Q^ 

(^1,^2) = l^^\\n\2^2 (39) 



and 

. [ei.6](i + ||CIP)-(6,C)[C,6]-[ei,C](C6) .... 



where r/i = (if7a(^i) and 772 = c^c7a(^2)- Basis vectors of the tangent space 
T^^(^)7q(M^") are just the images of basis vectors from the parameter space. 
We see that these scalar products are rather complicated in comparison with 
the analogous scalar products on the parameter space. 

Now let us consider a quadratic 'energy' function on the projective space 
P„(C) 

Em = (41) 

that can be locally presented as 

^(0 = ^(7.(01^^17.(0) (42) 

where H is some Hermitian operator on C""*"^, and suppose that H has 
as its non-degenerate eigenvector. -E(C) is a representation of the function 
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(HTj) within the chart Cq, = {Ua, fa, K^")- Standard Euchdean and symplectic 
inner products defined by Eqs. fl30l) - fl31l) exist on the parameter space M^". 
Since basis {ej,iej},j 7^ a is symplectic, matrix of the symplectic form fl3T!) 
is 

a<") = (_°^ (43) 

Without loss of generality in the remainder of this section we assume that 
a = and suppress index a in all forthcoming expressions. 
After realification the energy function E{x,y) becomes 



XiXjAij + E y^yjAij 2 ^ XiyjBij 

i,j=l i,j=l i,j=l 

2(1 + ||x||2 + \\y\ 



^(^' y) = on , ii^ii. , ll..l|2^ (44) 



where Eq is the eigenvalue of H corresponding to the eigenvector eo, and 
A = Re H, B = Im H are real and imaginary components of the operator 
H matrix with respect to the basis (1281) . A is symmetric and B is skew- 
symmetric. Partial first derivatives of E{x,y) are 

n n 

2xjE{x, y) -J2 ^ji^i + E Bjiyi 

y) = . , n^l. , n n.'^' (45a) 



dxj 1 + ||a;P + ||?/P 

n n 



"^VjEix, y) ~Y1 Ajiyi - E Bj 



Realified matrix of the second derivatives at the origin {x, y) = (0, 0) is 

A - Eoln -B 



Spectra of operator H — EqIu in the Hermitian space C" and the Hessian 
d^QQ^E in the Euclidean space M^" are identical. Indeed, it is easy to show 
that 

det [IK - iuhn] = I det [H - {Eq + cu)I„] 1^ (47) 

Eigenvalues are just the 'transition energies' uJio = E^ — Eq. Positive definite- 
ness of Hessian implies that the function fHTj) has its minimum at the point 
eo- 
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Now let us try to exploit the symplectic structure of the parameter space. 
In theory of symplectic manifolds any smooth real- valued function on a sym- 
plectic manifold is called a Hamiltonian function. In particular, both function 
(j^Tj) and its local representative fH2l) are the Hamiltonian ones. 

In general symplectic manifold is a pair (M, uj) where M is an 

even- dimensional differentiable manifold and uj is closed skew-symmetric 2- 
form on M that is, for each x G M the mapping '■ T^-M x T^^^M — > R is 
non-degenerate bilinear skew-symmetric, Ux varies smoothly in x, and du = 
{d is the exterior derivative of u). For example, the mapping ( — >■ u^^^ (see 
Eq. (l40!) ) endows the surface 7q,(]R^") with the symplectic structure. 

For any smooth function / : M — M its differential at point x G M is a 
covector d^f '■ T3.M — >■ M. The image of d^f with respect to the isomorphism 
(see Eq.(l5])) is a vector of the tangent space T^^M and the first order 
differential equation 

j^lit) = {d,it)f) (48) 

is called a Hamiltonian one. Critical points of / are just the singular points 
of the vector field x (d^f) and vice versa. 

In certain situations to study the stability of solution of Eq. (H8|l in a neigh- 
borhood of some its singular point it is sufficient to analyze the linearization 
of this differential equation (see, e.g. [IS]). 

For the energy function (jUj) its differential is a covector 

" d " d 

j=i ^ j=i y^ 

that can be transformed to the the symplectic gradient to give 

©J'W.,.)^) = -Ee,7^E(x,y) + X^ie,Ai^(x,y) (50) 

j=i yj j=i j 

(here the standard symplectic structure of the parameter space is used). 
The Hamiltonian equations in coordinate form are 

Xj= -^E{x,y) (51a) 
y^= ^E{x,y) (51b) 



where j = 1, 2, . . . , n. 
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Linearization of Hamiltonian equations in a neighborhood of its critical 
point means that Hamiltonian function is replaced by its quadratic approxi- 
mation. For E{x, y) in a neighborhood of the origin we have 



(52) 



J,i=i «,i=i «J=i 

The corresponding linearized Hamiltonian equations are 

n n 
Xj^ - E (A - EolnhVj - E ^ij^j 

n n 

Vj = E (A - EoIn)ijXj - E (53a) 
j=i j=i 

or, in a matrix form 

Characteristic roots of real matrix —Q'K are purely imaginary. Indeed, 

det + Xhn] - I det [H - {E^ + iA)I„] 1^ (55) 

and, consequently, the aforementioned roots are ±ia;jo. 

Realified parameter space (C")k ~ R^" can be again complexified to give 
complex vector space C (S>m IR^"" of complex dimension 2n with basis vectors 
{1 Cj, 1 ® icj} where 1 is the unit of C Hermitian structure on this space 
is introduced in the following way: for any a, a' £ C and any f , u' G R^" 

{a <^ v\a' ® v') = aa'{v,v') (56) 

where (•, •) is standard Euclidean scalar product on R^". Basis {1(8)6^, l<S>iej} 
is orthonormal with respect to this scalar product. 
Let us select in C (8)^ R^"" a new orthogonal basis 

fj = ^ [(1 ® e,) - i(l ® ie,)] , fj = ^ [(1 ® ej) + i(l ie,)] (57) 
The corresponding transformation matrix is 
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In this basis the Hamiltonian equations (151bp take the Schrodinger-type form 



C,= 2i4E(C,C) 

C,= -2i|-E(C,C) (59a) 



where (j = xj + iyj and where, by definition, 



d 1 f d d , , , 

- ^TT- (60a) 



dCj 2 \dxj dyj 
^ (60b) 



d(j 2 \dxj dyj 

On the complexified parameter space C <S)r M^" energy can be considered as 
a function of independent variables (, (. 

Linearized matrix Hamiltonian equation fl5^ in variables (,( becomes 




= -C-^n^KC ^ ^ (61) 
where 

C-'S)KC^-,(«-„^»'» -(H -%„.„)) e'^' 

Spectra of real Hamiltonian matrix —QO-C and complex Schrodinger-type ma- 
trix —C^^Q'KC are obviously identical. It is clear as well that the last matrix 
is diagonalizable over C. This means, in particular, that real non-symmetric 
matrix is also diagonalizable over C. Over the field M of real num- 

bers this matrix can be transformed into block-diagonal form with n skew- 
symmetric 2x2 blocks. Indeed, if Ui is an eigenvector of matrix —Q'K 
belonging to the eigenvalue icUjo then Ui is also the eigenvector of this ma- 
trix belonging to the eigenvalue — itUjo- It is easy to show that vectors ReUj 
and Imui constitute a basis of two-dimensional invariant subspace of matrix 
—Q'K in real parameter space M^". In this basis matrix —Q'K becomes a 
direct sum of 2 x 2 real matrices and differential equation fl54l) becomes a 
direct product of n equations 



Pi] = ( ^io\ fp, 



(63) 
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where Pi, qi are coordinates of real vector from the aforementioned two- 
dimensional subspace relative to the basis {Re^i, Imw,}. Thus, matrix —Q!K 
can not be diagonalized over M but can be transformed to the following simple 
('canonical') form 

















^20 








V 


















(64) 



Note that spectrum of matrix is purely imaginary without depen- 

dence on character of critical point of the energy function. Identification 
of the index of the critical point under consideration can be performed by 
means of analysis of inequalities |Aj| H-i^o > Eq. Phase curves of Hamiltonian 
systems behave differently for critical points of different index. In particular, 
the solution of linearized Hamiltonian system reasonably approximates the 
solution of the initial non-linear system if the Hessian of the Hamiltonian 
function is sign-definite (see, e.g., [T6]). 

It is pertinent to mention that by a certain abuse of notation we did not 
distinguish operator H defined on the space C""*"^ and its restriction on the 
subspace C" complementary to Ccq. 

Concluding this section we can state that for the projective spaces it 
is of no consequence what structure, Euclidean or symplectic, is used (if, 
of course, we are interested only in stability of energy critical points and 
excitation spectra but not in actual evolution). In the next sections it will 
be demonstrated that for submanifolds of projective spaces the situation is 
different: use of symplectic structure may give results essentially different 
from that obtained with the Euclidean structure. 



Hartree-Fock Manifolds 

We start with relevant assertions from the multilinear algebra. Their 
proof may be found, e.g., in [TTl HH |19]. Our presentation is close to that in 
[in]. The notion of the wedge product is supposed to be known. 
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Symbol JF^ will stand for 1-electron sector of the Fock space spanned by 
n = \N\ molecular spin orbitals (MSOs) {ipi} with indices from the MSO 
index set = {1, . . . ,n}. Vectors from JF^ = are called p— vectors 

by mathematicians and p— electron states by physicists, p-vector z is called 
decomposable if there exist vectors zi, . . . ,Zp e J-'j^ such that z = zi A . . . A 
Zp. In quantum chemistry decomposable p— vectors are called 'p— electron 
Slater determinants'. Interpretation of quantum chemical notions in terms 
of modern multilinear algebra may be found in [20] . 

Proposition 1. Vectors Zi, Z2, ■ ■ ■ , Zp from JF^ are linearly independent if 
and only if 

zi A . . . A 7^ 0. (65) 

Definition 5. For arbitrary p— vector z its annihilator is 

Ann{z) = {y e Tlf\z A y = 0} (66) 

It is clear that for any p— vector z its annihilator is a subspace of the 
one-electron vector space JF^. 

Proposition 2. Let z be p— vector with annihilator spanned by free vectors 
zi, Z2, . . . , Zg. Then there exists {p — q)— vector y such that 

z = Zi A Z2 A . . . A Zq Ay (67) 

In quantum chemistry annihilator of some p— electron state is called 'sub- 
space of inactive MSOs associated with this state'. 

Proposition 3. Let zi, Z2, . . . ,Zp and yi,y2, ■ ■ ■ ,yp be two free families of 
vectors from JF^. Then 

ZiA... Azp = XyiA... Ayp (A ^ 0) (68) 

if and only if p— planes (subspaces) generated by vectors zi, Z2, . . . , Zp and 
yi,y2, ■ ■ ■ ,yp are identical. 

Since Ann[z) = Ann{Xz) for any p— vector z and any A 7^ 0, it is possible 
to consider Ann as a mapping 

Ann : (K\{0})2;i A Z2 A . . . A Zp ^ [zi, Z2, . . . , Zp] (69) 

where 

p 

[zi,Z2,...,Zp] = '^Kzi (70) 
1=1 

is a K— linear hull of vectors zi, Z2, . . . , Zp. 
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Definition 6. The set of all p— dimensional subspaces (p— planes) of one- 
electron Fock space JF^ is called its Grassmann manifold and is denoted by 
the symbol Gp(jF^). 

It is easy to see that the mapping defined by Eq. fl69l) is actually a bijection 
and its inverse is an embedding of the Grassmann manifold Gp(jF^) into the 
the projective space P(JF^) of p— electron states. 

Definition 7. The set Ann~"'^(Gp(jF^)) is called the Hartree-Fock (HF) man- 
ifold of p— electron states. 

HF manifold can be characterized implicitly as the set of solutions of a 
system of homogeneous polynomial equations, that is as a projective algebraic 
variety (see, e.g., [21j). This characterization is based on the following simple 
statement (see, e.g., jlQj). 

Proposition 4. For any non-zero p— vector the dimension of its annihila- 
tor is less or equal p. Non-zero p— vector is decomposable if and only if its 
annihilator is of dimension p. 

Recasting the vector equation 

zAy = Q, zeJ'l, yeJ'], (71) 

in a coordinate form (with respect to some fixed one-electron basis) , we arrive 
at a homogeneous linear system of {J^^ scalar equations with respect to n 
unknowns 1/2, ... , and it is easy to see that p— vector z is decomposable 
if and only if all minors of order n — p + 1 of the matrix of this system are 
equal to zero. These conditions give us the required system of polynomial 
equations. 

Explicit characterization of HF manifolds in terms of local coordinates 
seems to be much more useful for applications. 

Let us fix some one-electron basis set {'i/'i,V'2, • • • ,'^n\ and consider K- 
linear hull of vectors 

n 

Zi = ^ipjZji, i = 1,2,..., p (72) 

To this linear hull the mapping Ann~^ puts into correspondence the line 
generated by decomposable p— vector 

z = ZiA...AZp= ^ V^i^ A . . . A (73) 

l<il<...<ip<n 
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where is the determinant of p x submatrix of n x p— matrix Z with 

row indices ii, . . . ,ip. 

Coordinate charts on Grassmann manifold may be introduced as follows. 
Let us suppose that Zi-^,,,i^ ^ 0. Then we can write 

Zi 

J2 



(2:1 ... , Zp) = {^l^i,, . . . , ^/^i^, Vji, . . . , ( 2 ) ( ''4) 



where Zi is non-degenerate p x p— submatrix of matrix Z with row indices 
ii, . . . ,ip and Z2 is its {n—p) xp— submatrix with complementary row indices 
ji, . . . ,jn-p- It is clear that vectors {zi . . . , Zp)Z]^^ generate the same p— plane 
[zi, . . . , Zp]. If another set 

(2/1 ■ ■ ■ , Vp) = {Ai, ^ip, ^ji, • • • , ( ) ^^^"^ 

of free vector generating this plane is chosen then necessarily Y2Yjf ^ = Z2Z^^. 
Thus, as a local parametrization it is possible to take the mapping 

lhi2...ip ■■ Z (^ii, . . . , ipip, ^ii, • • • , ^in-p) 2 ) ^ ^ ^2 A . . . A 2;^] (76) 

where Z G and 

n—p 

= + ^i^jiZnk, {k = l,2,...,p) (77) 

The domain of the corresponding chart is Ui-^i^,,,ip = (K^^""^-*). It is 

easy to see that the family of charts 

where ^Pij^i^...ip = forms an atlas of Grassmann manifold and that this 

manifold is K— analytic. Note as well that the chart with indices ii, . . . , is 
centered at the point (p— plane) [■i/'jj, . . . , ijji^. It is pertinent to mention that 
the described atlas of the Grassmann manifold depends on the choice of the 
MSO basis set and by a properly selected non-degenerate transformation of 
MSOs any point of this manifold can be placed at the center of the 'stan- 
dard' chart C12...P = (t^i2...p5 V^i2...p5 IK*''^""^)). For p = 1 Eq. fl76l) is identical to 
Eq.([2ni). 
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In addition to the aforementioned 'canonical' realization of the Grass- 
mann manifolds, there exist another realizations, of which we mention three 
most commonly used ones. 

(1) The set of all Hermitian idempotents (density operators) p over JF^ 
such that Tr p = p; 

(2) Quotient of the general linear group GL(^^) modulo its certain closed 
subgroup. 

(3) Quotient of the unitary group U(JF^) modulo its certain closed sub- 
group. 

Two last realization require additional explanations. There is a natural 
transitive action of the general linear group GL(^j^) on the set of all p— planes 
from J^lf. If TT G Gp{J-'jf) then its isotropy group G.„ = {g E GL{J^lf)\giT = n} 
is a closed subgroup of GL(J^^). If p first vectors of a chosen one-electron 
basis generate p— plane vr then matrix representation of transformation g e 
Gtt is of the form 

The quotient space GL(jF^)/G'7r can be endowed with K— analytic struc- 
ture consistent with the quotient topology [22]. For any p— plane vr the set 
GL(jF^)/G'^ of left cosets of the general linear group GL(jF^) relative to 
is a homogeneous space isomorphic to the homogeneous space Gp(jF^). 
Indeed, let us put 

0^ : gG^ ^ grc (80) 

We have 

(i) Onig'gGj,) = g'gn = g'^gn) = g'9T,{gG.„) for any non-degenerate trans- 
formation g' from GL(jF^); 

(ii) 9^{gG^) = 9^{g'G^) ^ gn = g'jr =^ gG^ = g'G^; 

(iii) For any n' G Gp(jF^) there exists g G GL(jF^) such that n' = gn ^ 
n' = e^{gG^). 

Property (i) means that the mapping is a morphism of homogeneous 
spaces. Properties (ii) and (iii) imply that this mapping is a bijection. 

Thus, for any fixed p— plane vr the Grassmann manifold Gp(JF^) can be 
identified with the space 

GL(^^)/G^= U gG^ (81) 
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of left cosets of the general linear group of one-electron Fock space relative 
to subgroup Gj,. 

The general linear group GL(jFj^) is a Lie group (that is a group and 
K— analytic manifold) and its tangent space = T/^(GL(JF^)) at the 

identity /„ is a Lie algebra of all one-electron linear transformations that is 

S(^^)=^^®(^^)* (82) 

Subgroup Gt, is also a Lie group and its Lie algebra Q{n) is a subspace of 
the vector space 0(JF^). And again, if p first vectors of a chosen one-electron 
basis generate p— plane vr then Lie algebra g(7r) can be identified with the 
algebra of matrices of the form of Eq. flT^ with elements from the number 
field under consideration. 

Let us consider a decomposition 

0(^^)=0(7r)©p (83) 

where subspace p is constituted by matrices of the form 

[It] (-) 

It can be proved (see, e.g., [22j) that there exists a neighborhood of zero in 
the parameter space where the mapping 

7^ : Z ^ exp (Z) G^ (85) 

with 

p n 

and 

6.;. = |^.)(^^| (87) 

is a K— analytic parametrization of a neighborhood of hi GL(jF^)/G7r. 
The family of mappings {77r}7r involves infinite number of members. Due 
to compactness of Gp(jF^) there exists a finite subfamily of this family 
that parametrizes this manifold. If vr is some p— plane spanned by MSOs 
ipi, ■ ■ ■ y'ipp then the mapping 

Tpa'^.p : Z - [exp(Z)^i] A ... A [exp(Z)^J (88) 
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may be used to parametrize representatives of HF lines in the projective 
space of p— electron states belonging to a neighborhood of p— plane vr. 

Both mappings fl76l) and fl88|) are K— analytic. But they have a certain 
drawback. Namely, even if the set of MSO corresponding to the origin is or- 
thonormal (with respect to the standard Euclidean or Hermitian scalar prod- 
uct on one-electron Fock space), the parametrized p— frames corresponding 
to non-zero values of parameters are not. This may be inconvenient both for 
evaluation of matrix elements and in the course of solution of optimization 
problem. It is easy to modify the definition of the mapping (!76|) to eliminate 
the aforementioned drawback: 

7?i...p : Z ^ (^1, . . . , V^„)g(Z) ^ [g(Z)V^i] A ... A [g(Z)^p] (89) 

where 




and W(Z) is px p upper triangle matrix performing Gram-Schmidt orthogo- 
nalization of vectors representing p— plane from t/i2...p, and where the initial 
one-electron basis is supposed to be orthonormal ( as has been already men- 
tioned, without loss of generality it is possible to consider p— planes from 
Ui2...p). For p = 1 the mapping (!89|) coincides with the mapping (1271) . For 
the case of the complex parameter space of dimension p{n — p) its realifi- 
cation leads to the real parameter space of dimension 2p{n — p) with basis 
{eu^^iCy^}. This basis is orthonormal with respect to the Euclidean scalar 
product (Z, Z') = Re [Tr Z'''Z'] and symplectic with respect to skew-symmetric 
scalar product [Z, Z'] = Im [Tr Z'I'Z']. 

To modify properly the definition of the parametrization fl88|) . it is nec- 
essary to consider the unitary subgroup U(JF^) of the general linear group. 
Unitary transformations possess the following important properties: (1) uni- 
tary group acts on the set of all p— planes transitively, and (2) if p— plane vr 
is invariant with respect to m G U(jFj^) then orthogonal complement vr^ is 
also invariant with respect to u. As a result, the isotropy group of arbitrary 
p— plane tt is a direct product Gt, = U(7r) x U(7r^). Lie algebra u(jF^) of 
the unitary group consists of skew-Hermitian matrices with, in general, com- 
plex elements but it is a vector space over the field M of real numbers (after 
multiplication by, say, the imaginary unit skew-Hermitian matrix becomes 
Hermitian). The parameter space for the exponential parametrization should 
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be taken as the space of all matrices of the form 



n p 



Z = ^ ^ [Zy^Cy^ - Z^f,e^^] (91) 

u=p+l fJ.=l 

and parametrization mapping 7}^ ^ is given by Eq. (l88|) but with parameter 
matrix fl9Tl) instead of the matrix fl86|) . 

Orthogonal with respect to the trace inner product basis in this real 
parameter space can be chosen as 

o-ufj. (^ufi (^fii/ (92a) 
Suf, = i (e^f, + (92b) 

In this basis the parameter matrix (19 ip takes the form 

n p 

Z= ^ 5^ [X.^a.;. + (93) 
1/—P+1 fj,=i 

where X, Y are real and imaginary components oi {n ~ p) x p parameter 
matrix Z. 

Symplectic form on real even-dimensional parameter space with basis 
Wufi, Sufj.} is defined with the aid of coordinate functionals rfX^^ : Z Xy^ 
and d/Yi/^ '. Z — > as 

n p 
u=p+l fi=l 

If 7i2...p(Z) is a local parametrization of representatives of HF states (built 
with orthonormal MSOs) then 2p{n — p) vectors 

||_.(Z)^g.......|L. ...... (-.) 

^^(Z) = X:z,A...A^A...Az, (96b) 

Ufl j^-j^ Ufl 

form a basis of the tangent space to the HF manifold at the point 712.. .p(Z). 
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We start with the energy functional parametrized by 7i2...p(Z): 

^(Z) = ^(7il.p(Z)|i/|7il.p(Z)) (96) 

where if is a Hermitian operator acting on the p— electron sector JF^ of the 
Fock space. We suppose that this functional reaches its minimum at the 
origin X = Y = of the parameter space and the corresponding single- 
determinant wave function is $o = V^i A . . . A ipp. 

When using exponential parametrization 7}^2 we confine ourselves 

to the quadratic approximation of energy functional and, consequently, to 
the linearized version of the Hamiltonian equations. 

Local expansion of 7}^2 p(Z) in a neighborhood of the origin is 



iTr (ZtZ) 



\,i A, A' i<i' 



% + *'^A. + E E •^'/^A.^A'.' + ■ ■ ■ (97) 



From this expansion it readily follows that basis of the tangent space to 
the HF manifold at the origin is constituted by 2p{n — p) 'single excited 
determinants': 



ax.,. dY,. 



(0) = $;, ^^(0) = i$;: (98) 



Energy differential at the origin is 



n p 

dio,o)E = - E J2[^K\H\'^o) + {<^o\H\<!>;)]dX,, 

v=p+l tJ.=l 

+ ^ E E + ($o|i^|i$P] dY^, (99) 

and it is easy to see that stationary conditions (i(o,o)-E'(ay^) = d{Qfi)E{sy^) = 
are equivalent to the well-known in quantum chemistry Brillouin conditions 
{^'^\H\^o) = 0. Note that we have not specified yet the concrete Hermitian 
operator involved in Eq. fl96l) and, consequently, even the classic form of the 
Brillouin theorem is of a rather general nature. 

At this stage it is convenient to introduce the following matrices: 

A.^,.V = a-S.u')il-5^^.'){%\H\<^'';;,) (100b) 
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Note that H*^^^ is Hermitian whereas A is symmetric. 

Differential of the quadratic part of energy function (jHS]) at a point (X, Y) 
of the parameter space is 

rf(x,Y)^^^^ = {dX dY) :H (101) 

where 

_ /Re (H^is + A) - Eolp(n-p) -Im (H^is + a) \ 

^ " V (RCis - A) Re (RCis _ A) - i5;oIp(n-p)y' ^ 

is the matrix of the energy function Hessian calculated at the origin of the 
parameter space with respect to the basis (!92|) . Here Eq = {^q\H\^q) = 
2E{0). 

The next step is to use the symplectic form (1941) to get the linearized 
Hamiltonian equations analogous to Eqs. flS^ for the case of the projective 
manifolds. We have 




where 



n = ( ° ^ (104) 

\ — Ip(n-p) / 



is the matrix of the symplectic form fl9^ 

Complexification of the parameter space gives 2p(r;,— p)-dimensional com- 
plex space where the matrix Hamiltonian equation (I103p takes the form 



H^is _ Eolpin-p) A \ f Z 

-A -{m^ - Eolp(^r.-p)) J 



(105) 



At this stage we have three matrices: real symmetric matrix of the second 
order derivatives !K, real Hamiltonian matrix —Q'K, and complex Schrodinger 
matrix —C~^Q'KC where 

C = - (yi^-P^ (log) 

2 y ilp(n— p) ilp(?i— p) J 

First two matrices are relative to the basis fl92|) of real parameter space and 
the third matrix is relative to C— transformed basis {1 (S> o-jy^, 1 (S> s^,^} of 
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the complexified parameter space. Using isomorphism 'i(o,o)7i2...p5 "^-e can 
identify the parameter space with the tangent space to the HF manifold at 
the point ^(0, 0). With such an identification matrices "K and —Q'K can 
be considered relative to the basis of single excited determinants fl98|) whereas 
Schrodinger matrix becomes relative to the orthogonal basis 

T;: = ^[(l®<fP + i(l®i$a (107b) 

of the complexified tangent space. 

Matrix "K of the second order derivatives is orthogonally diagonalized 
over R and has real spectrum consisting of 2p{n — p), in general case dif- 
ferent, eigenvalues. Hamiltonian matrix QJi obviously satisfies the equality 
[Q'KYQ + QIQJ-C) = which means that it is an infinitesimal-symplectic ma- 
trix (element of Lie algebra of the symplectic group Sp(M^'^)). It is easy to 
show that spectrum of such (non-degenerate) matrix satisfies the following 
restrictions: (i) if A is its eigenvalue then necessarily —A is also its eigen- 
value; (ii) each complex eigenvalue appears in pair with its conjugate. And of 
course, Hamiltonian and Schrodinger matrices have identical spectra. These 
matrices, however, are neither Hermitian nor skew-Hermitian. They are not, 
in general, orthogonally diagonalized. In contrast to the case of the projective 
space, spectrum of these matrices is not necessarily purely imaginary. 

Till now p— electron Hermitian operator H was not specified and all equa- 
tions obtained are valid for arbitrary choice of this operator. For electronic 
Hamiltonian all matrix elements involved, say, in Eq. fllOSp . may be easily 
calculated with the aid of the standard technique to give 

H^^^V = [Eo + eu-e^]6f,^>6:,^' + {ufi'Wfiu') (108a) 
A.^,.'^' = ififi'lWi^') (108b) 

where (ij II A;/) = {ij\kl) — {ij\lk). Substitution of these expressions in Eq. (11051) 
immediately gives the standard TDHF equations (see, e.g., p]). 

In analogy with classical mechanics, matrix A in Eqs. fll03p - fll05p may be 
interpreted as a constraining matrix and its appearance is of the same nature 
as appearance of constraining force in Newton equations. The case A = 
corresponds to CIS method which reduces to diagonalization of operator H 
projection on the subspace of p— electron states spanned by vectors {^Z}. 
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As we have already seen, choice of Euchdean or symplectic metric on the 
projective space leads, in essence, to equivalent theories. For HF manifold, 
however, it is not the case. And in our opinion there is no a priroi crite- 
rion, what metric is preferable. Strictly speaking, it is necessary to compare 
behavior of exact excitations energies, TD excitation energies, and Hessian 
eigenvalues as functions of (complex) parameters A^^y^/ to make reliable 
conclusion. In Appendix A analysis of behavior of excitation energies is 
performed for a simple model case. 

Now let us turn to Gram-Schmidt parametrization (l89l) . Using Eqs. 
(B.15)-(B.16) from Appendix B, it is easy to get the following general ex- 
pressions for partial derivatives: 



d 



dX, 
d 



-7if..p(Z) = J][W(Z)],,$^^(Z)-Re[g(Z)Wt(Z)l $(Z) (109a 



i/fi 



7l2...p 



(Z) = 5^[W(Z)]^,i<l>^^(Z) + ilm [g(Z)W^(Z)j __^ i$(Z) (109b) 



where 



$(Z) = ZiAz2A...AZp = 7i2 ..p(Z) 
$J(Z) = zi A . .. A Zj^i Aipy A Zj+i ... AZp 



;iio) 
'111) 



Energy derivatives are easily calculated at arbitrary point Z = X + iY. In re- 
alified parameter space ]R2p("~p) symplectic scalar product may be introduced 
as the imaginary part of the Hermitian trace product. After transformation 
of energy differential to symplectic gradient it is possible to write down the 
exact Hamiltonian equations: 



X 



UfJ, 



Y 



2Im 



2Re 



^iy,,(z)(<i.(z)|/j|<i.^^(z)) + (g(z)w^(z))^^E(z 



p 



J2w,,{Z)mZ)\H\^'<iZ)) - (g(Z)W^(Z))_ E(Z) 



(112a) 
(112b) 



where E{Z) = |($(Z)|i7|$(Z)). 
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The corresponding Schrodinger-type equations are 



Zy^ — — 2i 



5^W^^,(Z)($(Z)|i/|$J(Z)) - (g(Z)Wt(Z)) E{Z) 

X;iy,,(Z)($(Z)|i7|$J(Z)) - (g(Z)Wt(Z))_J?(Z) 
.i=i 



(113a) 
(113b) 



Remind once again that these equations are exact (not hnearized) ones. Their 
hnearization will not lead to new equations, because it does not depend on 
the choice of local parametrization (see, e.g., [T6]). 

The most common in quantum theory approach uses Grassmann mani- 
fold without its prior embedding into p-electron projective space. This ap- 
proach is applicable to general energy functionals that can be defined in 
terms of 1-density idempotent operators. Since exponential parametrization 
of Grassmann manifolds can be found (in explicit or implicit form) in almost 
all pubhcations concerning TD theories (see, e.g.,|2l [3]), we confine ourselves 
to analysis of Gram-Schmidt parametrization. Instead of the mapping ( l89i) 
we have 



TT 



GS 
12. ..p 



Z ^ ^ \zk){zk\ 



(114) 



fc=i 



where Zk = g{Z)ipk and g(Z) is given by Eq.(|90l). To simplify notations, the 
explicit indication on dependence of matrices g and W on parameters X, Y 
will be omitted. For the same reason we suppress subscript 12. . .p (chart 
multiindex) and superscript GS of the parametrization mapping flll4p . 

Expressions for partial derivatives of the mapping 7r(X, Y) are easily de- 
rived from Eqs.(B.6)-(B.9) of Appendix B: 



d 



9X 



-vr 



Ufl 



(X, Y) = J2 [Wfik\i^u){zk\ + W,k\zk){i^u\] 



k=l 

p 



;il5a) 



d 



k,l=l 
P 



-IT 



(X, Y) = i^ [W^km{zk\ - WMiM 



k=l 
P 



i ^ [Wf^ig^k - W^kQui] \zk){ 



(115b) 



kJ=l 
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In particular, real tangent space to the Grassmann manifold at the origin is 
spanned by the vectors 

d 

7r(0, 0) = e^f, + e^jy (116a) 



-|-7r(0,0)= i(e,^-e^,) (116b) 

where e^^ = 0:^,1, fi = l,...,p, and u = p + 1, . . . ,n. 

The second order derivatives of the parametrization mapping at the origin 
are also easily calculated to give 

7r(0, 0) = 6^^' {tyyi + tyiu) " ^ vv> (c^^' + 6^'^) (l 17a) 



7r(0, 0) = -iS^^'{euu' - eu'u) - ^5uu'{e^^: - e^/^,) (117b) 



dXufidYu'f,' 



vr(0,0)= 7r(0,0) (117c) 



Let us suppose that orthonormal MSO basis is fixed and each 1-electron 
operator is identified with its matrix. In particular, Grassmann manifold 
Gp(JF^) can be identified with n x n— matrices p satisfying the following 
restrictions: 

p'f = p, Trp =p, = p (118) 

A function E{p), smooth with respect to real variables aij = Re pij, j3ij = 
Im Pij and, in general, complex- valued, will be referred to as 'the energy 
function', or just 'the energy'. Within the chart under consideration (with 
multiindex 12 ... p) local representative of the energy function restriction to 
the Grassmann manifold is £(X, Y) = EorrQi, Y). It is reasonable to suppose 
that for physically relevant energy functions the imaginary component of 
£(X, Y) vanishes (note that, even if this condition is fulfilled, the imaginary 
part of E{p) should not necessarily be equal to zero for arbitrary complex 
matrix p). To avoid cumbersome expressions, we suppose that the energy is 
a holomorphic function of complex variables pij. In this case realification of 
the energy domain is not required. 

Using Eqs. flllSp . it is easy to calculate partial derivatives of £(X, Y) at 
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arbitrary point (X, Y) of the parameter space: 



(9X, 



■(X,Y) 



(X,Y) = i 



(In - ggO ^ gWt 
(In - gg^) ^^gWt 



— 1 



Vfl 



Wgt^ (l„-ggt) 

Wg^^' (in-ggO 



(119a) 



(119b) 



where ^ is, in general complex, matrix of partial derivatives calculated 
at the point 7r(X, Y). 

The Schrodinger-type evolution equations on the complexified parameter 
space are 



' UfJ, 



2i 
-2i 



Wgtf^ (l,-ggt) 
"(l„-ggt) 0gW^ 



(120a) 
(120b) 



To linearize these equations in a neighborhood of the origin it is neces- 
sary to calculate the first and the second derivatives of £(X, Y) at the point 
(0,0). It can be easily done with the aid of Eqs. flll6l) - flll7p . The linearized 
Schrodinger-type equations are 

d^E , dE r dE 

+ O^^i- 0,j,ji- 



'Z'up, — '^ 



E 



dp, 



a^-E 



dp^,udp^>yi 



;i21a) 



U'fJ.' 



d'^E , dE 



dpuf^dpi^'u' dp 

vv' 

d^E 



dE 



dp, 



MAI 



Zu'fl' 



(121b) 



Derivatives on the right-hand side of these equations are taken at the point 

Classic HF theory supplies us with an example of simple energy function 
that can be defined in two ways: 



4'f(p) = 2Trp[h + F(p)] 



122) 
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or 

41(p) = ^Trpt[h + F(p)] (123) 

where 

F^j{p) = K^ + Y,Pikm\3l) (124) 

k,l 

These functions coincide on the subspace of Hermitian matrices p but differ- 
ent as functions on the space of all complex matrices. In particular, E^{p) 
is holomorphic as a function of complex variables whereas E\^^[p) is not. 

It is easy to show that Fock matrix F(p) is Hermitian for any Hermi- 
tian p and that the energy (and, consequently, E^) restriction to the 
Grassmann manifold is a real- valued function. 

We confine ourselves to the holomorphic case. Simple calculation give 

^^(p) = F,,(p) (125a) 



dpij 



^HF 



(p) = {jl\\ik) (125b) 



dpijdpki 

where, within the chosen chart, p(X, Y) = 7r(X, Y). 

The Schrodinger-type evolution equations on the complexified parameter 
space, corresponding to the energy function E^{p), are is readily obtained 
from Eq.([T20D 



Z,, = 2i [WgtF(p)(l„-ggt)]^^ (126a) 
Z,, = -2i [WgtF(p) (I„ - ggt)]^^ (126b) 

Using formulas fll21b )- fll21b ). it is easy to derive the linearized version 
of these equations: 

Zufi = i^ [{up'Wpu') + 5uu'5f,f,'{eu - £f,)] ^i.v 

+i^{uu'\\pp')Z^>^> (127a) 

-i'y^^{pp'\\vv')Zu'^i (127b) 
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where it is assumed that the origin of the parameter space is placed at a 
stationary point of £|^p(X, Y) and the canonical HF MSOs corresponding to 
this point are selected. 

The energy E^^p (non-holomorphic case) should be treated either as a 
function of real variables aij and Pij or complex variables pij and pij. 

Conclusion 

Physical TD theories, put properly in the framework of modern differen- 
tial geometry, may become a general and powerful tool for investigation of 
many electron systems. As soon as general scheme of derivation of evolu- 
tion equations in the case of, say, relatively simple complex projective spaces 
is elaborated, the same scheme with minor technical modifications can be 
applied for any projective algebraic manifold. Of course, existence of ad- 
ditional mathematical structures on the manifold under consideration may 
result in plenty of different equivalent realization of the same geometrical 
object and each realization may carry its own unique feature. Selection of 
relevant realization depends, of course, on the concrete task. In present pa- 
per we considered complex Grassmann manifold and its realizations as (1) a 
submanifold of many electron projective space which we called HF manifold, 
and (2) the set of idempotent 1-density operators. In our opinion, both real- 
izations supplement each other. Study of another realizations exploiting the 
fact that Grassmann manifolds are also homogeneous spaces of certain Lie 
groups may be found in [3]- [5]. 

Comparison of linearized Schrodinger evolution equations on many elec- 
tron projective space and on its HF submanifold readily reveals the appear- 
ance of a constraining matrix A in the HF case. This matrix includes ex- 
plicitly double excitations from the HF state and characterizes in a certain 
sense the curvature of the HF manifold. Perturbation analysis of eigenvalues 
(TD transition energies) of the matrix of linearized Schrodinger evolution 
equation on HF manifold leads to the conclusion that these eigenvalues are 
concave functions of matrix elements of the constraining matrix whereas the 
exact transition energies are usually convex functions of the same parame- 
ters. This means that, if a given CIS transition energy is greater than the 
exact one then, at least in a neighborhood of A = 0, TD transition energy 
should be more close to the exact value than the CIS one. This feature of TD 
transition energies is not so clearly seen if realization of Grassmann manifold 
as the set of 1-density idempotents is used. 
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Appendix A. Simple Example of 
Excitation Energies Behavior 

Let us consider simple case corresponding to the projection of p— electron 
operator H on the subspacc spanned by determinants 

where $o = A ^"2 A . . . A ■^p. General form of relevant operator H matrix 
in this case is 

/ Eo A \ 
H(A) = H^is A {A.l) 
V At At B / 

where Ai^^y^/ is a matrix of free complex parameters, A and B are arbitrary 
fixed complex and Hermitian matrices, respectively. 
Writing matrix H(A) in the form 

H,(A) = H(0) +£W(A) (^.2) 

where 

= EE [^-.--vl^o)«'| + A.^,.vK')(*o|] {A.3) 

and treating W(A) as perturbation, we come to the following second order 
expression for the exact transition energies: 



u;,o(A) = u;,o(0) + E M (1 + ^Ji) (^-4) 

where are the eigenvectors of matrix H(0)(^o — ^o)- It is clear that 
a;io(A) are convex quadratic functions. 
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Matrix of the Schrodinger-type equation (11051) (the imaginary unit pref- 
actor is omitted) with respect to the basis (I1U7I) of the complexified tangent 
space may be written in the form (A. 2) with 

and (non-Hermitian) perturbation 

W{^) = 2Q{K) - 2g(A)t (A6) 

where 

vv' flfl' 

Prefactor 2 in Eq.(A.6) appears due to somewhat non-standard normahzation 
of basis functions (^07^: (T;;|Tp = (T;;|Tp = i. 

Zero-order matrix H(0) has real spectrum ±u;™(0) symmetric with re- 
spect to zero. If is eigenvector of block (1,1) of H(0) corresponding to 
the eigenvalue a;™(0) (CIS transition energy) then is the eigenvector of 
block (2,2) of this matrix with the eigenvalue —uj™{0). We confine ourselves 
to the perturbation analysis of non-negative part of matrix H(A) spectrum. 
The second order expression for transition energies in this case is 



, ,TD/, \ _ , ,TD 



7^ +^.0 



From (A. 8) it readily follows that, at least in some neighborhood of A = 0, 
TD transition energies are concave functions of parameters Ay^y^i. As a 
result, if for a fixed index i graphics of functions uj.iq{A) and C(j™(A) have 
crossing points, then at this points TD excitation energy takes on the exact 
value. If ujJ^PiGi) > cc;io(0) then there exists a certain neighborhood of the 
point A = where TD transition energy u;™(A) more close to the exact 
values than the CIS energy. Probably in this sense one should interpret 
frequently occurring in literature statement that TDHF approach is superior 
in many aspects than the CIS one (see, e.g. [5]). 

To visualize the exact dependence of transition energies on single param- 
eter K, let us take the concrete complex Hermitian 6x6 matrix relative to 
the basis {$o,$?,$t$i,$^,$?^} 

15 k- 

-11.8433 -0.31349-0.47024i 0.04167+0.1528i 0.2+i 
Vi(,^\ I -0.3135+0.47024i -9.9623 -0.54167-0.18061 0.3 1 ( A Ci\ 

= I -7 0.1-i I l^-yj 

0.04167-0.1528i -0.5417+0.1806i -8.1944 0.2-i 

R 0.2-i 0.3 0.1+i 0.2+i -f 
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and use graphical tools of the Mathematica package [23| . In FigH] depen- 
dence of the exact transition energies, TD transition energies and Hessian 
eigenvalues on real parameter n is displayed. 




Figure 1: Transition energies as functions of parameter k 

Exact transition energies are convex and TD transition energies are con- 
cave functions of k in certain neighborhoods of the origin for all four ex- 
cited states. The aforementioned intersection points are essentially differ- 
ent for different excited states. For example, ujJPiK) > ojfk^^ and 

"^acPl^) i^so^'^*- In Figl2] the second order approximations of TD 

transition energies are shown together with the exact and TD transition 
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energies. It is seen that at least in the concrete case under consideration the 




Figure 2: Exact transition energies uj'?^'^^, TD transition energies cu™, and 
second order TD transition energies ujfa'^^ as functions of parameter k 

second order perturbation correction (A. 8) gives excellent approximation of 
TD transition energies for all states. 

In some cases use of perturbation theory to evaluate excitation energies 
may be computationally more efficient than the direct diagonalization of ei- 
ther non-symmetric Hamiltonian real matrix or non-Hermitian Schrodinger- 
type complex matrix. However, to make reliable conclusion about applica- 
bility of perturbation approach, heavy numerical testing is required. 
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Appendix B. Analytic Expressions for 
Derivatives of Gram- Schmidt 
Parametrization Mapping 



Matrix 

g(Z) =(^\^ W(Z) {B.l) 

involved in Eg. (189!) satisfies the orthogonality condition 

[g(Z)]tg(Z) = l, (5.2) 

To simplify the notations, we agree to omit, when possible, explicit indication 
on dependence of matrices W and g on Z = X+iY. Differentiation of Eq.(B.2) 
gives 

1^ [Wt]-^ + W-^ = -Wt Jt^g - gt J,,W (5.3a) 



9Wt 



[WT ' + W-i^ = iWtjt^g - igtj.^W (5.36) 



where basis matrices Jjy^ and iJ,^^ are defined by 

= Jut. (BAa) 

- iJ.^ (5.46) 



d 


(z 






( 


'I, 




Csl 



Following Garton let us introduce the matrices 

P., = W-^^ (5.5a) 

Q.M = W-^^ (5.56) 

where u = p + 1, . . . ,n and /x = 1, . . . ,p. Since W is upper triangle, matrices 
and Qufj. are also upper triangle and their matrix elements are 

{0 if i > j 

-i[Wtjt^.g]^.^-l[g%,W]^^ if^ = j (5.6a) 
-[Wtjt^g],,-[g%.W]^^. if^<j 
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{0 if ^ > j 

i[WtjU],,-i[g%.W].. iii<j 

At the origin P,^(0) = Q,^(0) = 0. 
Prom Eqs.(B.5)-(B.6) it follows that 



dW 

= WP,^ {B.7a) 



dW 



gy- = WQ'^'^ (^-^^^ 
and 

= iJ.^W + gQ,^ {BM) 

Now it is easy to calculate the derivatives of vectors Zj — g{'Zi)ipj with 
respect to parameters Xj^^ and Y^fj^: 



3 



dzj 



92; 

^^"^^ k=i 



At the origin we have 



dW dW 



3z ■ 



Oz ■ 

^^^^ 
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